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Abstract 

In nature self-organized systems as flock of birds, school of fishes or herd of sheeps have to deal with 
the presence of external agents such as predators or leaders which modify their internal dynamic. 
Such situations take into account a large number of individuals with their own social behavior which 
interact with a few number of other individuals acting as external point source forces. Starting 
from the microscopic description we derive the kinetic model through a mean-field limit and finally 
the macroscopic system through a suitable hydrodynamic limit. 
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1. Introduction 

The aim of this paper is to present different level of descriptions for the dynamic of a large group 
of agents influenced by a small number of external agent. In a biological context this corresponds to 
the behavior of a flock or a school of fishes attacked by one or more predators, or the movement of 
a heard of sheep guided by a sheepdog. Recently such dynamics have been studied also in robotic 
research, where engineers tried to control the action of a school of fishes introducing a fishbot 
recognized as leader, see [1]. 

From the modeling viewpoint this turns out in considering a microscopic dynamic described by 
classical flocking models like Cuker-Smale and D ' Orsogna-Bertozzi et al. (see 1(3, 12]) interacting 
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with a set of few individuals characterized in an analogous way as in [9]. Moreover, motivated by 
the works of [H], we endowed the classical dynamic of interaction both with a metric as well as a 
topological interaction rule. 

Following the approach in [8| we start from the microscopic dynamic, given by a ODEs system, 
and we derive two other different level of description: the mesoscopic (or kinetic) level through a 
mean-field limit and the macroscopic level through a suitable hydrodynamic limit. At difference 
to the first order macroscopic models proposed in j9j] here we obtain second order models for the 
corresponding continuum dynamic. 

Finally we report some numerical examples for the solution of the mean field model in a series 
of test cases. The simulations have been performed using the fast algorithm recently presented 
in 1|. Rigorous mathematical results concerning the asymptotic limits just described can be found 
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2. Microscopic model 



We are interested in the study of a dynamical system composed of N individuals and N p 
external agents with the following general structure 

±i = Vi i — 1, . . . , N 

1 1 Np 

' il = JV* ^ F ( X i> V U X 3> V j) + JfZ-s^iXiiVuPk) (1) 
J6^(i.) P fc=l 

Ph = <Ph(t,p,Ap N (ph)) h = l,...,N p 

where (x, v), = (xi,Vi) lives in M. 2d , d > 1, i = 1, . . . , iV and (p)^ = ph € M. nd , with n = 1,2, 
h = l,...,N p , and i\T p < JV. 

Function i 7 describes the interactions inside the swarm, and -F p depicts the interaction with 
each external agent ph- According to the three zone model [3], F can be decomposed in 

F(xi,Vi,Xj,Vj) = H(xi,Xj)(vj - Vi) + A(xi,Xj) + R(xi,Xj) + (2) 

characterize the alignment term, ^4 the attraction, R the repulsion and S 1 represent a seZ/ 
propulsion-friction term. The same decomposition holds for F p if p^ = i.e. n = 2. 

In a first order model, n = 1, similar interaction can be considered, bu for the first term we should 
speak of consensus dynamic rather then alignment. 

Moreover we endow the model with a topological rule of interaction. Each agent will interact 
only with a fix number of agents of their specie, 

E^(xj) = {the N* closest neighbors respect to i} . 

A topological interaction is motivated by some recent studies, see 0, Qj]], suggesting that in a 
flock each individual modifies its position according to few individuals directly surrounding it, no 
matter how close or how far away those agents are. If N* = N each agent interacts with all the 
others and the topological interaction coincides with the global metric interaction [H, 0] • Functions 
ifh : [0, +oo) x (M. nd ) N p x R — > K, describe the evolution in time of each external individual and 
they depend on the discrete density p N defined as empirical measure 

1 N 

i=l 

According to 0] we define A as a convolutional operator Ap N (t) = p N (t) * rj, where rj is a smooth 
kernel with compact support. 

2.1. Classical swarming models 

The classical swarming models take in account a global interaction between the agents, that 
corresponds in our case to choose N* = N neighbors. 
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Cucker-Smale model. Describes an alignment dynamic, see and [3], 

H{xi 1 x J ){v j - Vi) = H(\xi - xj\)(vj - Vi), (3) 

where H(\xi — Xj\) is a function that measures the strength of the interaction between individuals 
i and j , and depends on the mutual distance, under the assumption that closer individuals have 
more influence than the far distance ones, as depicted by 

H(r) ' 



(l + r 2 )7' 

where 7 > discriminates the behavior of the solution, see 

0,3 and @. 



D'Orsogna, Bertozzi et al. , considers a self-propelling, attraction and repulsion dynamic, see 121 ] 
A(xi, Xj) + R(xi,Xj) + S(vi)vi = -V Xi W(\xj —Xi\) + (a — /3\vi\ 2 )vi (4) 

where W ■ M. d — > R is a given potential modeling the short-range repulsion and long-range 
attraction and a, (3 are positive parameters. A possible description is given by the following power 
law 

W(r) = -, 

a 

where a > b > are positive parameters. 

Both the models take in account symmetric interaction between agents, which correspond to 
the conservation of the first momentum. This assumption isn't really realistic if we want to model 
group of animals, to improve the models other features have been introduced as a perception 
cone [l|, 0] or the concept of relative distance 14 ] . In these cases we lose the interaction symmetry 
and consequently the conservation of the first momentum, note that also the topological interaction 
breaks the symmetry. 



3. Kinetic model 

In order to retrieve a mesoscopic description of the system (pQ) we proceed through a mean-field 
limit, as already done in 0, ; S|- Since the limit is done only for the first set of equations, which 
describe the evolution of the swarm, we obtain an hybrid model made by one kinetic equation and 
the ODEs system governing the external agents. 

The basic idea of the mean-field limit is to derive, through a weak argument, a single evolu- 
tionary equation for / , the empirical measures defined as 

1 - 

f N (x, v,t) = -J2 S(x - Xi (t))S(v - Vi(t)). 
i=i 

A second step is to show that the limit f N — > f for N — > 00 holds true, then rigorous derivation 
of the limiting kinetic equation can be proved. A well-posedness theory for this derivation has been 
developed in [7[ for a general set of swarming models. 
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We remark that in the mean-field limit the analogous of the topological set, of interaction 

is described by a characteristic function on the ball B(x,R*), with center x and radius R* such 
that 

R* = min I R 



R>0,p* = / f(x,v)dxdv } . 

jR d JB(x,R) J 

The model at mesoscopic level reads 

dtf + v- V x f = -V„ • (£*(x, v)f) - V, • (£ p (x, v)f). 

(5) 

ph = cph(t,p,Ap(p h )) h = l,...,N p 

where p(x,t) = f Rd f(x,v,t)dv. 

Remark 1. If we consider the decompositions H|) in the case of Cucker-Smale and D'Orsogna- 
Bertozzi et al. models and = (x?, v?), we have a more particular formulation for £* and £ p , 
this way 

£*(x,v)=— / H(\x - y\)(w - v)f(y,w)dydw+ 

P JR d JB(x,R*) 
1 

and 
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(6) 



V x W{\x-y\)p{y)dy + S(v)v 

P JB(x,R*) 



1 N p N p 

£ p & V ) = W H HP d x ~ - V )- W Y1 V.W*(|s - x\\) (7) 

p k=i p k=i 

4. Hydrodynamic model 

We are interested also in the macroscopic description of the system, from a numerical point of 
view this corresponds to reduce the dimensionality of the problem in such a way that simulations 
become affordable. As done in [1, [l3| we define 

, = *.,«) = //*., ^=*fl=/* T:=T( M ) = /k-»l^, 

in order to obtain a system of equations which describe the evolution of these quantities, we 
integrate the kinetic equation in ([5]) against dv and vdv . 

According to 0] we impose the closure the momentum assuming that the fluctuations are 
negligible, i.e., that the temperature T(x,t) = 0, and the velocity distribution is monokinetic: 
f(x,v,t) = p(x,t)5(v — u(x,t)). The previous assumptions lead to the following hydrodynamic 
system 

d t p + div x (pu) = 0, 

d t (pu) + V x ■ (pu ® u) = J 7 *(x,u)p(x,t) + J 7p (x,u)p(x,t) 1 (8) 
t Ph = (ph(t,p,Ap(p h )) h=l,...,N p 
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Figure 1: Shepherd dogs. With parameters a = 2.5, b = 0.1, 7 = 0.45 for models (O and (j4|. 



In the particular case of Cucker-Smale and D'Orsogna Bertozzi et al. model and ph = (x^,i>^) the 
second term reads 



F*(x,u)=— [ H(\x - y\)(u(y,t) - u(x,t))p(y,t)dy+ 

P JB(x,R*) 



B(x,R*) 

V x W{\x - y\)p(y, t)dy + S(u)u 



1 

P* Jl3(x,R*) 



(9) 



«) = W E ^(1* - ^DK - *)) - E V,^ P (k - 41) ( 10 ) 
JV P fc=i JV P fc=i 

5. Numerical examples 

We show here two numerical examples taken from [d| for two different biological dynamics. In 
order to solve the kinetic part of the model, the simulations are performed using the Asymptotic 
Binary Interaction algorithms proposed in We recall that the algorithm is based on a stochastic 
routine which use N s sample particles. The overall cost is 0(N S ). 
The dynamics considered for the swarm in both the cases are 

F(xi,Vi,Xj,Vj) = H(\xi - Xj\)(vj - Vi) + V Xi W(\xi - xj\), 

and only a repulsion dynamic respect to ph in F p , given by the W p (r) = —r c /c with c > 0, and 
where r = \x% — Ph\- 

Shepherd dogs. The simulation represents the evolution of a swarm controlled by two leaders pi,P2 £ 
R 2n , n = 1, interacting with the swarm in the following way 

<Ph(t,p,8 h ) = V p u h \ = ; V p = 300, r p = 5, h = l,2. 
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= ^6 ( max {°' r P - l x l 2 }) 2 ' s h ■= Ap(x p h ) = (p * x Vr])(x p h ). 
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Figure 2: Swarm attacked by a predator. With parameters a = 4, b — 2, 7 = 0.45 for models ((3]) and 



Swarm attacked by predator. We consider the evolution of a swarm which undergoes the action of 
a predator. The predator is modeled by the evolution of p = (x p , v p ) £ M 2n , n = 2 and its evolution 
is lead by the following potential 
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